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Abstract. We investigate a classical mass-action model for diseases transmitted by arbovirus. 
Most traditional studies have focused on global stability issues and oscillatory behaviour. Here, 
we take a different view and we study the dynamics of the model, under the hypothesis that the 
typical timescale for the host and the vector are distinct. In this way, two asymptotic dynamics 
naturally appear: the fast host dynamics and the fast vector dynamics. The former yields, at 
leading order, a SIR model for the hosts, but with a modified incidence rate. Thus, the vector 
disappears from the model, and the dynamics is similar to a directly contagious disease. The 
latter yields a SI model for the vectors, with the hosts disappearing from the model. Numerical 
results show the performance of the approximation, and a rigorous proof validates the reduced 
models. 



1. Introduction 

1.1. Background. Vector-borne diseases are different from direct contagious ones since there is 
indirect transmission from host to the vector and vice- versa. From the point of view of mass-action 
epidemiological modelling, perhaps the most natural and simplest model is the coupling of SIR 
model for the host and SI model for the vector that was first developed by [1, 3 . It is given 
schematically in Figure [l] 
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Figure 1. Compartmental description of the arbovirus model by [fllS]. 
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Parameter 


Meaning 


and 


Number of hosts and vectors 


vl and 


birth rate for hosts and vectors 


(31 and (3: 


probabihty of a host being infected by a vector and vice-versa. 


7* 


recovering rate 


m* 


Number of alternative blood sources 



Table 1. Description of parameters meaning in the compartmental model de- 
picted in Figure ITl 



The most straightforward translation from the compartmental model shown in Figure [T] is given 
by the following system: 



(1) 



Q* /I* Q* (^v b C* 7"* 

Oy — ^ fJ^yOy ^^^^^ Dyl^ 

r* _ f^yb o* r* _ 



Here starred variables indicate dimensional quantities. The variables /S, / and R have the usual 
epidemiological meaning, with the subscript indicating if they refer to the hosts or to the vectors. 
For the other parameters, we refer to Table [l] 

System ([T]) has been comprehensively studied by [4^ , and extensively used for studies of Dengue 
as described, for instance, in [5]. Most theoretical studies have focused in the global stability 
features of System ([T]). In particular a proof of global stability using the theory of competitive 
systems can be found in [4 . More recently, some partial Lyapunov proofs were also available in [2] 
and |7] . Also a proof using a Lyapunov approach for the disease free equilibrium, and a persistence 
argument for the endemic equilibrium can be found in |8j. 

In this study, we shall take a different route and first show that there are two natural and 
independent timescales in System [l] namely the host timescale for the disease evolution and the 
correspondingly one for the vector. We then examine the possibility that these two scales are 
asymptotically separated and find that, for fast vector dynamics. System ([T]) can be approximated 
by a 2-D SIR model, with rational incidence rates. For fast host dynamics, the simplification is 
even more dramatic, and System ([T]) can be approximated by a 1-D SI model, again with a rational 
incidence rate. The derived approximations are uniform for all positive times. Incidentally, some 
results on global stability are also presented. 

1.2. Outline. In Section[2j we describe the basic model studied and review some of its properties. 
A global analysis using only a Lyapunov approach is also presented. In Section [3j we study the 
Fast Vector Dynamics limit, i.e., 

cr, /i^ ~ e~ , 

with all the remaining nondimensional parameters 0(1). This yields a reduced SIR system with 
modified nonlinear incidence rates. The dual limit, the Fast Host Dynamics^ i.e. 

with all the remaining nondimensional parameters 0(1) yields a reduced SI system that is some- 
what less interesting. The latter limit is presented briefly in AppendixjAj Conclusions are presented 
in Section m 

2. Preliminaries 

2.1. Scalings and adimensionalisation. We adimensionalise System [l] by letting 
{SI II Rl) = N*^{Sh, h, Rh) and {S;, ID = /„). 
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Also 

e = (A*)"^t, 

where A* is, for the moment, an arbitrary time-scale. We immediately obtain the new System 

Sh = - Sh) - SShlv 
ih = ^Shlv - il^h + l)h 
Rh = ih - l^hRh 

Sv = — Sy) — CrSylh 

ly dSylfi l^yly 

where 



/j.*A*{N;1 + m*) A*(7V* + m*) 

Also 

7* 1^1 . 1^1 

Typically, A* will be chosen such that at least one of 7, /i^, are order one. We shall not make 
any a priori choice, but indicate appropriate choices when studying particular regimes. 

As showed in [4 , the dynamics is restricted to the nonnegative orthant, and it is conservative: 

Lemma 1. Let R>q denote the nonnegative orthant ofR^. Then a solution of System with 
an initial condition in M>q^ stays in M>q. Moreover, the sums -\- -\- and S* + /* are 
preserved by the evolution. 

Because of Lemma [T] if the initial conditions for the host fractions add to one, with the same 
being true for the vector fractions initial conditions, then this will be preserved by the evolution. 
Therefore, we can work with the simplified, but equivalent, model below: 

r X = fih{i-x)-sxz 

(2) { ^ = SXZ-{^h^j)Y 

[ Z = a{l- Z)Y - flyZ 

2.2. Global stability analysis. System Q has the following two equilibrium points 

(1) The disease free equilibrium: X* = 1, F* = Z* = 0. 

(2) The endemic equilibrium: 

^* ^ il^h + 7)m^ + l^hcr _ lih [^cr - (/i/, + 7)/i^] _ /i/, [8(7 - {/j^h + j)lJ^v] 



The endemic equilibrium can be conveniently written as 

^* ^ ^ 1 + RqDq ^ Rq -'^ Z* = D 



where 



i^O 1 + ^0 ' /i/, + 7 ^0 + i^O^O ' ^l + i?o^o 

Ro = — 7 ^ — r and I^^o = 

or as 

= J^l±RoDo * ^ — ^(l-X*), Z* =I)o(l-X*) 
Ro l^Do ' /i/. + 7 

Theorem 1. Let Rq be defined as above. Then for Rq < 1 the disease-free equilibrium is globally 
asymptotic stable, while for Rq > 1 the endemic equilibrium is globally asymptotic stable. 

Proof. Suppose Rq < I and consider the following Lyapunov function 

v(x,r,z) = x-iogX + y' + ^z. 

fjjy 
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Then 



On substituting 



v = xl 1-- ) + z 



V = -fih ^^ - ifih + 7)(1 - Ro)Y - RoZY, 



which is negative for < i^o ^ 1? and X, F, Z in the interior of the positive octant of M^, other 
than = (1,0,0). 

For Ro > 1, let 

V{X,Y,Z)=X-X*log—+Y-Y*log— + (z-Z*loi 



Z* 



Then we have 



1 + -X 



X* 



■ - (/^/^ + 7)^ — (1 - z)- 



3-X* -X 



X* 
X~ 



(1 - 



crJX* 



(/i/^ + 7) (i^o^*(l - Z*) - 1) F 



cr^X*Z* Y dXZy 



fly Z 



Y 



Since 



we are left with 



3 - X* - X 



X* 



aSX* cr(5X*Z* Y 6XZY* 
Y — 



Writing 

and on using that 



fly fly Z 

(X*)2 X*(l-X*) 



Y 



X 



X + > 2X* 

X 



we find 



This can be rewritten as 



aSX 



YZ-i^hil-X*) 

fly ys. 



X fly fly Z Y 

X* fih(j{l - X*) Y 5fik{l - X*) XZ 



By the Arithmetic- Geometric inequality, we have 
Hence, 



- X*) y 5^i,,{l-X*)XZ ^ ^ ,o. ..*^l/3 



< -SUhil - X*) {RoX*f/\ 



V < 3/ih(l - X*) (l - (i?o^*)'/^) - "^^YZ. 



Since 



^0-^ = , , ^ — > 1, smce Rq > 1, 



1 + i^o 
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we find that V < in tiie interior. □ 

3. The fast vector dynamics 

If we assume that the vector population typical times are much sorter than the typical times 
for the host, we might expect that the vector population will be nearly in equilibrium, and hence 
that Z = 0. In this case, we formally obtain from Q 



(3) 



Z = 



Equation (|3| can be recast as a SIR system with a modified incidence rate. While the derivation 
above is heuristic, we shall see that it can obtained from a consistent multiscale asymptotic 
expansion. Moreover, the formal expansion can be rigorously justified. In what follows, we shall 
write 

a = d-e~^ and /i^ = Jlye~^ . 

Then, we are interested in solve 

r X = ^h{i-x)-sxz 

(4) { ^ = SXZ-{^h^-f)Y 

[ eZ = a{l- Z)Y - fLyZ 

subject to the initial condition 

X(0)=Xo, Y{0) = yo and Z(0) = Zq. 
3.1. Asymptotic expansion. Let 

er = t. 

Then, we seek an expansion of the form 

X = XO(t) + +0(e) 
Y = Y\t)^0{e) and 
Z = ZO(t) + ZO(r) + 0(e), 

where 

lim Z°{t) = 0. 

The leading order equations are given by the following differential-algebraic system: 

X^ =fih{l-X'')-5X''z'' 
Y,''=SX^Z^-{^h^j)Y^ 
0=a{l-Z^)Y^-jiiyZ^. 



This yields 

and hence, we obtain the system: 
(5) 



with initial condition X^(0) = Xq and Y^{0) = Iq- To solve for Z^, we notice that 

X^{t) = X\er) = X\0) + erX^iO) + 0(6^), 

with similar expansions for Y^(t) and Z^(t). 
Therefore, we find that Z^ satisfies 

ZO = -(<7Fo + M.)^°, 
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I.e 

Z\t) =(zo- J^' ) e-(^^o+^^)^ 

V (TYo lily J 

3.2. Global stability analysis of the asymptotic system. First, we have the following basic 
result 

Proposition 1. Let 

S = {x -\- y < 1^ x^y > 0}. 

Then S is invariant by the flow of In particular, the corresponding solutions are global in 
time. 

Proof. Since = is an invariant set, solution with ^ > will remain so. Also, when = 0, 
the flow points inside S. Thus x>0. Finally, {X^ + Y^)t < 0. □ 

System ([5| has two equilibrium points in the positive quadrant: 

(1) X* = 1, = 0. 
(2) 

Remark 1. Notice that X^ = X* and = F*. Thus the equilibrium of ^ correspond 
to the projections in the XY plane of the equilibria of Notice also that 



We then have 
Theorem 2. Let 

Ro 



Then Rq = Rq, ct'i^d for Ro < I the disease free equilibrium is globally asymptotically stable. For 
Rq > 1, the endemic equilibrium is globally asymptotically stable. 

Proof. Let 

F(Fo) = 

Then F is a Dulac function for system ([5| in proper compact subset of 5, since 

dxo IM^ - Xo)F{Yo) - SXo] + dy, [SXo - {im + 7)(^>o + M.) = -^hFiY^) - 6 - + 7) < 0. 

Thus, the system cannot have a closed orbit in the interior of S. For i^o ^ 1, the only equilibrium 
in 5 is (1,0). Thus by all orbits must converge to this equilibrium point. 
The linearisation of ([5| is 




For the disease free equilibrium, the eigenvalues of Jacobian are — /i/^ and {jj^h +7)(i?o — !)• Thus, 
the disease free equihbrium is a locally asymptotically stable node for Rq < 1^ and a saddle for 
Rq > 1. The unique orbit that approaches the disease free it easily shown to be the intersection 
of S with Yq = 0. Thus, all the other orbits must approach the endemic equilibrium. 

□ 
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3.3. Asymptotic convergence. The asymptotic expansions derived in |3.1| can be shown to be 
indeed asymptotic. The ideas used here are similar to the ones used to formaUse Kinetic Menton's 
theory — cf. [6 and references therein for instance. In particular, we have the following: 

Theorems. Let W{t) = {X{t),Y{t), Z{t)) andW^{t,r) = {X^{t),Y^{t),Z^{t)^Z^{r)). Denote 
the uniform norm in [0, oo) by ||.||oo- Then, for sufficient small e > 0, there exists a constant 
C > 0, independent of e, such that 

WW -WHoo<Ce. 



The proof of Theorem [3] is divided in several lemmas. 
We write system ([5| in a more concise form as 

W = J'(W,Z), 

with W(t) = {X{t)^Y{t)y and and Q being the appropriate entries of the right hand side of 
Q. We shah write 

W = W V eW° + eQ and Z = Z V ^ V eZ, 

with 

Notice that since is bounded, we need only to prove that (Q, Z) exist, are bounded, and are 
unique. In this case, we can then take C = ||(Q, ^)||oo- First, we observe that 

WO + W = J'(W^ Z^ + Z^) + K{t) and + ^'^ = g{W^^ Z^ + Z^) + 

where 

K{t) = SZ\t/e){X^{t) - Xo) and L{t) = a Z^ {t / e) {Y^ {t) - Yo). 

In particular, because of the fast decay of Z^, and because 1/(0) = K{0) =0, it follows that there 
exists a constant C > 0, such that 

POO POO 

/ K{t)dt, / L{t)dt<Ce. 
Jo Jo 

Since J" and Q are quadratic, we write: 

T{W, Z) =7"(W°, Z° + + eDw-^(W°, Z° + Z°)Q + eDz J"(W°, Z° + Z°)Z+ 
+ eHXY (^J 

a(W, =g(W°, Z° + + eDw6^(W°, Z° + Z°)Q + eDza(W°, Z° + Z°)Z+ 
where Q= 

Let T{t, s) be the fundamental solution to the linearised system 
I.. /Dw-^(WO, + Z^) DzHW°, Z^ + /Q 

Then direct integration yields 
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Lemma 2. The functions (Q, ^) satisfy the following integral equation: 

1 



/t(m)(^ 



(7) - / T{t,s) [^y'> ] As- - j T{t,s) I 1-1/ I ds. 




Moreover, the last term is bounded uniformly in e. 

We also have the fohowing large time behaviour result for the linearised system (|6|: 
Lemma 3. Let {Q{t)^Z{t)) a solution to (|6|. Then 

lim (Q(t),Z(t)) = 0. 

In particular, the solutions to (|6| are bounded uniformly in time any given e. Moreover, they are 
also uniformly bounded in e < 1, for all t > 0. 

Proof. For notation convenience, let us write ([6| as 

Fix e > and {Xq, Yo,Zo). Prom Theorem [2) we know that 

lim (W°(t), Z°{t) + Z°{t/e) = {W*,Z*), 

t^OO 

where (W*, Z*) will be the globally asymptotic stable equilibrium given by Theorem [l] Therefore, 
there exists T > 0, such that t > T implies that A is negative-definite. Since ^ 0, as e ^ 0. 
We can choose T such this holds for all < e < 1. 

Because of the continuity of the solution with respect to the initial conditions, we have that T 
is a continuous function of the initial conditions. Since these lie on a compact set, we can pick T 
such that A is negative definite for all < e < 1 and for all {Xq^Yq^ Zq)^ such that Xq^Yq > 0, 
-^0 + ^0 ^ 1 and < Zq < 1. .But then, for any such (Xq, Yq^Zq) and < e < 1, we have that 

lim T(t,T) = 0. 

t^oo 

Therefore, for any initial condition {Qo^RoY^ we shall have 



□ 



hm T(t, 0) = hm T(t, T)T(T, 0) = 0. 



Proof of Theorem^ First, we observe that the nonlinear term in ^ is locally Lipschitz, hence a 
standard fixed point yields existence and uniqueness for < t < to , for some, possibly small, t^. 

From the decaying of and from Lemma [ij we conclude that the two last terms of ^ are 
uniformly bounded in time and e. 

Moreover, Lemma [s] implies that the first two terms on the right hand side of ^ are also 
uniformly bounded in time and in e, for sufficiently small e. Thus, we conclude that the same also 
holds true for (Q, i?)^ Therefore, the solution to ^ is globally defined in time, and it is bounded 
uniformly in e, if the latter is sufficiently small. 

□ 
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3.4. Numerical results. We now present some numerical results in Figure [2] There we compare 
the full model with the asymptotic model. As expected, the approximation of X{t) by X^{t) and 
of Y(t) by Y^{t) are indeed uniform for all time, while the approximation of Z{t) by Z^(t) fail 
to be uniform in an initial layer. Notice also that such nonuniform behaviour is suppressed by 
including the corrector term in the initial layer. 

Remark 2. We observe that the structure of the asymptotic solution shows that the solutions must 
be confined, except for an initial layer, on what is usually termed a slow manifold. This behaviour 
is indicated in Figure^ 



4. Concluding remarks 

Diseases that are vector-borne have a number of features that distinguish them of contagious 
ones. Typically, timescales for the dynamics of the host and vectors are not within the same order, 
since mosquitos, for instance, that are a prevalent vector for such diseases have a very fast life cycle 
compared to humans say. With this in mind, we investigated the dynamical consequences of having 
host and vector dynamics with distinct timescales in the classical arbovirus model introduced by 
[Ul^. The natural regimes to look in this model are the fast vector dynamics (FVD) and fast host 
dynamics (FHD). While the former seems to be the most natural choice, we take the view that 
there might be scenarios where the latter may be observed. 

By means of a formal multiscale asymptotic analysis, we study both regimes. For the FVD, we 
find the leading order dynamics yields a SIR model for the host, with a modified incidence rate. 
Thus, the vector is removed from the model being present only parametrically as a function of 
the host infected fraction. Such a relationship, appart from its mathematical interest, might also 
be useful in verifying if field data is conformant, within the model, with the regime hypothesis. 
Additionally, the FHD regime yields an even more dramatic reduction with an SI model for 
the vectors, again with a modified incidence rate. Numerical results presented show that the 
approximation is indeed uniformly asymptotic in time. 

An interesting feature of the studied regimes is that they do not imply any condition on Rq^ 
and hence are compatible with a variety of disease developments from the point of view of global 
dynamics. Indeed, for both reduced models, the equilibria are preserved by the asymptotic approx- 
imation, and the global stability dynamics is consistent with the global stability dynamics of the 
full model. Finally, we have confirmed rigorously the asymptotic character of the approximation 
up to the derived order. 

The results presented here suggest that other epidemiological models of vector-borne diseases 
might be studied in a similar fashion. Hopefully, this procedure would provide model reductions 
allowing a better understanding and, perhaps, even a classification scheme for different arboviruses 
diseases based on the relative timescales of infection, the host and the vector dynamics. 



Appendix A. The fast host dynamics 

Analogous to the previous case, we now assume that the dynamics of the host is much faster 
than the dynamics of the vector. In this case, one might expect that host populations is nearly in 
equilibrium. Hence, we should have X ^ and F ~ 0, i.e, we should have the system 

[ Z = (7(1- Z)Y - iiyZ 

Equation ([8| can be seen as SI system for the vector, with a modified incidence rate. 
In order to justify System ([8|, we assume that 

5 = Je""*^, jih = Uh^~^ and 7 = 76""*^. 
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Figure 2. Comparison between the full and asymptotic models. Here Rq ~ 1.58, 
Do 0.75, e 0.01. Also, j^h 0.3, S 0.4, 7 0.35, a 0.5 and j2y 0.2. 
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Figure 4. Slow manifold dynamics 
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Figure 5. Comparison between the full and asymptotic models for the suscepti- 
ble and infected host fractions. Parameters are the same parameters, except for 
jj^h = 0.005 and 7 = 0.4. In this case, we have that Rq = 2.47, Dq = 0.0125. 



Thus, we are interested in solve 

r eX = ilh(l-X)-5XZ 

(9) < eY = dXZ-{fih^j)Y 

[ Z = Z)Y - ii^Z 

subject to the initial condition 

X(0)=Xo, F(0) = ro and Z(0) = Zq. 

In what follows, we formally derive the leading order asymptotic expansion and provide a global 
analysis together with some numerical results. The proof of asymptotness is very similar to the 
FVD regime, and hence it is omitted. 

A.l. Asymptotic expansion. As before, we let er = t. The asymptotic expansion now take the 
following form: 

X{t) = X\t) ^ X\t) ^ Q{e), 
Y{t) = Y\t)^Y\T)^0{e), 
Z{t) = Z\t)^0{e). 

Here, we shall also have 

lim (XO(r),f^(r) = 0. 
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At leading order, we have 

= j2h{l - X°) - 5X°Z° 
= 5X0^0 - iPh + l)Y° 
Z° = a{l - Z°)Y° - ^,Z° 

We can solve for and obtaining 

SZ + ph + 1 SZ° + flh 



Thus the last equation becomes 



Also, we write 

Z^{t) = Z°{eT) = Z°{0) + erZO(O) + 0{e^T^), 
hence, since Z°{0) = Zq, we obtain 

1° = flhX° - SX°Zo 

= SX°Zo = (7 + Ph)Y° 

The solution can be written as 

It is straightforward to verify that the eigenvalues of A always have negative real part, and hence 
that 

lim (X°(r),y°(T)) =0. 



A. 2. Global stability analysis. The equilibria are Zq = and Zq = Z*. Since 



dZn 



^oDq— 1 - RqDq ^ J^\2^0 

^0 + Zo [Do + ZoY 



At Zo = 0, its value is i^o — 1- So the origin is globally asymptotically stable for Ro < I. 
At Zo = Z* its value is 

1 + i^o^o 



i^o(l + ^o) 



(1 - Ro) 



Hence Z* is globally asymptotically stable, if Ro > I. 
When i?o = 1, we have 



Zo 



RoDoi^-l 
Do^Zo 



DoZo 



Do 



(Zo + 1) <0, Zo >0. 



Hence, Zo = is also globally asymptotic stable when Ro = 1. 



A. 3. Numerical results. The results for the components are qualitatively similar to the fast 
vector dynamics, and hence are omitted. Nevertheless, the reduction of the dynamics to the slow 
manifold is more dramatic in this case as shown in figure [6] 
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Figure 6. Slow manifold dynamics for the fast host regime. 
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